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We present first-principles calculations of the structural phase behavior of the [1:1] 
PbTiOa/PbZrOs superlattice and the PbTiOs and PbZrC>3 parent compounds as a function of 
in-plane epitaxial strain. A symmetry analysis is used to identify the phases and clarify how they 
arise from an interplay between different kinds of structural distortions, including out-of-plane and 
in-plane polar modes, rotation of oxygen octahedra around out-of-plane or in-plane axes, and an 
anti-polar mode. Symmetry-allowed intermode couplings are identified and used to elucidate the na- 
ture of the observed phase transitions. For the minimum-period [1:1] PbTi03/PbZr03 superlattice, 
we identify a sequence of three transitions that occur as the in-plane lattice constant is increased. 
All four of the phases involve substantial oxygen octahedral rotations, and an antipolar distortion 
is important in the high-tensile-strain phase. Inclusion of these distortions is found to be crucial for 
an accurate determination of the phase boundaries. 

PACS numbers: 77.80.bn,77.55.hj, 77.55. Px,77.84.Cg,68.65.Cd 



I. INTRODUCTION 

Ferroelectric materials have long been the subject of in- 
tense study both because of their fundamental scientific 
interest and because their properties are promising for de- 
vice applications. Experimental and theoretical studies 
of the relationships between their symmetries, structures 
and properties are essential for the development of im- 
proved materials. One of the most widely used and thor- 
oughly studied materials is PbZr x Tii_ x 03 (PZT), which 
in bulk shows an enhanced piezoelectric effect near the 
rhombohcdral-tetragonal morphotropic phase boundary 
(MPB) around x = 0.5. Nohcda et ali explained this 
anomalous response in connection with their discovery of 
a narrow wedge of monoclinic phase, bridging between 
tetragonal and rhombohedral phases at the MPB. 

Modern growth techniques such as pulsed laser deposi- 
tion and molecular-beam epitaxy now allow for the out- 
of-equilibrium synthesis of atomic-scale perovskite super- 
lattices of arbitrary layer sequences, and control of the 
in-plane lattice constant via coherent epitaxy on chosen 
perovskite substrates^ As the behavior of a bulk solid 
solution and a superlattice having the same overall com- 
position can in principle be quite different, these develop- 
ments offer an exciting opportunity for the development 
of improved new materials. However, given the enormous 
number of possible structures that can be constructed in 
this way, it is clear that theory has an important role to 
play in guiding the search. 

Not surprisingly, then, there has been a simultaneous 
explosive growth in the application of first-principles the- 
oretical methods to the study of perovskite superlattice 
structures^ specifically to the prediction of structure, po- 
larization and polarization-related properties such as the 



dielectric and piezoelectric response. Initially, the fo- 
cus was on the electrostatic interaction of the different 
layers, specifically on how a ferroelectric instability in 
one constituent could induce polarization in the other 
constituent^ 5 - maintaining a uniform displacement field 
throughout.— The systems considered included superlat- 
tices of PbTi0 3 (PTO) and PbZr0 3 (PZO), specifically 
the [1:1] PTO/PZO superlattice with varying epitaxial 
straini 7 - However, it has long been recognized^ that oxy- 
gen octahedron rotations play an important role in de- 
termining structure, polarization and related properties 
in ferroelectric perovskites. In first-principles work, such 
antiferrodistortive instabilities were identified in a num- 
ber of ferroelectric perovskites^ and it was shown that 
octahedral rotations in PZT play an important role in 
the morphotropic phase boundar y 11 ' 12 and possibly in 
high-pressure phases of PTO Most recently, combina- 
tions of rotations at the interfaces were found in first- 
principles investigations to generate improper ferroelec- 
tricity in PTO/SrTi03 superlattices^ and in layered per- 
ovskite compounds! 15 ' 16 An overview of the effect of rota- 
tional distortions on perovskites can be found in Ref. ^ 

In the present work, we apply first-principles DFT 
methods to determine the ground-state phases of the [1:1] 
PTO/PZO superlattice as a function of epitaxial strain, 
as well as those of pure PTO and PZO for reference. We 
go beyond previous work 7 - by including the octahedral 
rotations and tilts (both of which shall henceforth be re- 
ferred to as "rotations" ) which can be accommodated in 
a \[2 x \[2 x 2 supcrcell; the types of distortions identified 
as playing an important role are shown in Fig. [TJ We find 
that octahedral rotations are ubiquitous throughout the 
phase diagram of the PTO/PZO superlattice, and have 
an significant influence on the polarization and the loca- 
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FIG. 1: Sketch of distortions of type (a) Pool, (b) Pno, (c) Rooi, (d) Rno and (e) Aj 10 in the PTO/PZO superlattice. 



tions of the phase boundaries. In particular, they help 
stabilize a monoclinic phase over a very broad range of 
epitaxial strain. We relate these results to the behavior 
of the PTO and PZO constituents and analyze the re- 
sults in terms of the symmetries of the various phases; 
our approach is analogous to that for A-site superlattices 
reported in Ref. [l6|. Our results clarify the role played by 
the coupling between these degrees of freedom in deter- 
mining the structure and polarization of the superlattice 
structure under consideration. 

This paper is organized as follows. In Sec. [TT| we de- 
scribe the methods used to apply DFT to bulk PTO and 
PZO and to the PTO/PZO superlattice. The results of 
these calculations are presented is Sec. IIIIl and a discus- 
sion is given in Sec. IIV1 We summarize and conclude in 
Sec. El 

II. METHODS 
A. Computational details 

The effect of rotations on the PTO/PZO superlat- 
tice was studied using density-functional theory (DFT) 
calculations within the local density-approximation 
(LDA)i 18 i 19 The calculations were performed using the 
Vienna Ab initio Simulation Package (VASP)3i2i The 
Brillouin zone was sampled with a 4x4x4 Monkhorst- 
Pack gridi22 The projector augmented wave method 
(PAW)2i^ was used with an 800 eV plane-wave cutoff. 
The polarization for all systems was calculated using the 
Berry-phase method^ as implemented in VASP. 

For comparison with the PTO/PZO superlattice, we 
performed analogous studies of PTO and PZO with the 
same epitaxial constraints. For all three systems, the 
structure was first initialized by imposing selected po- 
lar and/or rotational distortions, not limited to those 
shown in Fig. [T] and the system was then relaxed sub- 
ject to the constraint of fixed in-plane lattice parame- 
ters, corresponding to the epitaxial boundary condition. 
While the symmetry of some of the phases allowed for 
a tilting of the c axis away from the normal to the a-b 
plane, tests described below in Sec. Ill Dl show that this 
effect is small, and the calculations reported here were 
all carried out with the c axis kept parallel to z. The 



in-plane lattice parameters were varied over the range 
3.80 A < a < 4.10 A for PTO, 3.90 A < a < 4.20 A for 
PZO and 3.90 A < a < 4.25 A for the PTO/PZO super- 
lattices. 

In order to consider the effects of rotations of oxy- 
gen octahedra, the calculations were carried out on a 
20-atom \/2 x \/2 x 2 unit cell. This supercell allows 
distortions with wavevectors at the T point (0,0,0), X 
point (0, 0,7r/a), R point (7r/a, 7t/a, tt/cl) and M point 
(%/a, n/a,0). In particular, it allows R point and M point 
oxygen octahedron rotational distortions which are rep- 
resented most generally by aTb~c~ and a~b~c + in the 
Glazer notation^ 

In all cases considered, we found that structures with 
Pioo > had higher energies than those with Pno> 0. 
Therefore in the following we limit our consideration to 
structures with nonzero Pool and/or Pno as shown in 
Figs. QJa-b). Furthermore, when comparing PTO and 
PZO phases incorporating i?-point (out-of-phase) rota- 
tions to those with M-point (in-phase) rotations, we al- 
ways found the i?-point patterns to be lower in energy. 
Similarly, with one exception that will be noted in due 
course, we always found the rotations of the ZrOe and 
TiOg octahedra to have opposite signs in the PTO/PZO 
superlattices. In what follows, therefore, out-of-phase ro- 
tations of the two layers stacked along z are to be under- 
stood unless otherwise specified. Rotations about [001] 
and [110] axes will be denoted by Rooi and Rno respec- 
tively, as in Figs.QJc-d). The 20-atom supercell also ac- 
commodates antipolar distortions of the type shown in 
Fig. HJc), denoted as Aj 10 . Structural energies for the 
PTO, PZO and PTO/PZO systems are given relative 
to the energy of the corresponding ideal high-symmetry 
P4/mmm structure (i.e., with no rotations or polariza- 
tion) at a reference lattice constant of 3.99 A. 



B. Small distortions 

Typically we started our structural minimizations from 
an initial geometry of relatively low symmetry, corre- 
sponding to a mixture of several kinds of distortion, so as 
to be conservative in our assumptions about what kind 
of structural distortions would actually occur. For ex- 
ample, in the PTO calculations to be described below, 
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we started from a system with the Pooij Piio> Rooi and 
Rno distortions of Figs. Q^a-d) present. After structural 
relaxation, some of these components (especially the ro- 
tational ones) would become very small but not disappear 
completely. To determine whether these small distortions 
were physically meaningful, or just a result of incomplete 
convergence or related numerical difficulties, the distor- 
tions were manually turned off (i.e., the structure was 
symmetrized) and then the total energy was calculated as 
the distortions were added back in in small increments. 
If the resulting energy curve showed a minimum away 
from zero distortion amplitude, the distortion was taken 
to be physically meaningful. This analysis proved useful, 
in particular, in clarifying the phase diagram for PTO in 

sec.nnsi 

C. Location of phase boundaries 

As we shall show later, the case of the PTO/PZO su- 
perlattice presents a related difficulty. Here, we find two 
second-order transitions at which one or more order pa- 
rameters go to zero in an apparent square-root singu- 
larity. In order to get a more precise location for the 
critical epitaxial lattice constant at which this occurs, 
we have carried out the following procedure. We de- 
scribe this method for the case of the transition from 
the lower-symmetry Plal space group to the higher- 
symmetry Plan space group with increasing epitaxial 
lattice constant in the tensile-strain region of PTO/PZO, 
but similar considerations apply to the transition from 
Plal to PAmm in the compressivc-strain region. 

First, we chose a lattice constant near the transition, 
but far enough away that the distortion leading to the 
low-symmetry phase was well developed. In our exam- 
ple, we chose a=4.125A, in the Plal phase but near 
the transition. We symmetrized this structure using the 
symmetry operations of the higher-symmetry phase, con- 
structed the eigenvector in displacement space pointing 
from the symmetrized to the low-symmetry structure, 
and used this to define a symmetry-breaking mode ampli- 
tude. Next, for each lattice constant near the transition 
(here for a=4.100 to 4.250A in increments of 0.025A), 
we started from the relaxed high-symmetry structure and 
added distortions corresponding to this eigenvector with 
several amplitudes (and no further relaxation), fitting to 
an expression of the form 

AE(u) = ^ku 2 + ^au 4 (1) 

to the energy vs. mode amplitude u (relative to u=0) 
as computed from the first-principles calculations. Not 
surprisingly, we found that k is strongly dependent on in- 
plane lattice constant a, while a is not. We then plotted 
k vs. a (typically obtaining a nearly linear behavior) and 
identified the critical lattice constant a cr it as the one at 
which this curve crosses through zero. As we shall see 
later, we found the transition to occur at a cr i t =4.17lA 



for the transition in this example. As we shall also see, 
using the same analysis on the compressive train side 
we found the transition from Plal to PAbm to occur at 
a crit =3.890A. 



D. Treatment of c-axis tilting 

As mentioned above, some of the space groups consid- 
ered in this work, specifically the monoclinic ones, allow 
by symmetry for the c axis to tilt away from the [001] 
direction, typically by developing a small [110] compo- 
nent. To study the effect of such a shear strain, especially 
on the location of second-order phase boundaries as de- 
scribed in the previous subsection, we again made use of 
the distortion pattern leading from the higher-symmetry 
tetragonal or orthorhombic state to the lower-symmetry 
monoclinic one as obtained in Sec. Ill CI Again denot- 
ing the amplitude of this normalized distortion by it, and 
starting from the higher-symmetry structure at the lat- 
tice constant of interest, we carried out calculations of 
changes in total energy and stress as a function of u and 
shear strain (axis tilt) n. We then parametrized the re- 
sults to the first-principles calculations using a general- 
ized version of Eq. (JXJ) taking the form 

AE(u, 77) = ^ku 2 + iau 4 + /3un + \pi? ■ (2) 

Again k was the only parameter that was found to 
depend strongly on a. Setting dE(u,n)/dn = and 
dE{u 1 rj)/du = within the model, we can estimate both 
the shift of the phase boundary and the amount of c-axis 
tilt. For the phase boundary when the PTO/PZO su- 
perlattice system goes from the Plal space group to the 
P2an space group, for example, this resulted in a shift 
of the critical in-plane lattice constant by only 0.0016 A, 
to 4.173 A, leading to small shear strains near the phase 
boundary (~0.2%). Similar results were found for other 
phase boundaries between monoclinic and tetragonal or 
orthorhombic phases. 

The impact of the tilting of the c axis appears to be 
less profound in our study than in that of Ref. [D, in which 
octahedral rotations were not included. Indeed, we find 
it to be small enough that it will not change any of the 
results significantly, and so we have not included it in the 
results to be presented in the next section. 



III. RESULTS 

A. Parent phases 

The results of the DFT calculations on PTO are given 
in Fig. [2j For in-plane lattice parameters less than about 
a = 3.90 A, we find the stable structure to be the tetrag- 
onal PAmm structure (space group 99) in which only 
Pooi is present (ferroelectric "c phase," 5-atom primitive 
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FIG. 2: (Color online) Calculated properties of PTO as a 
function of epitaxially constrained in-plane lattice constant, 
(a) Minimum-energy structure: octahedral rotations (black) 
of type Rno (*); polarizations (red) of type Pooi (°) and 
Pno (x). (b) Total energy (black squares) and c/a ratio (red 
triangles). 



FIG. 3: (Color online) Calculated properties of PZO as a 
function of epitaxially constrained in-plane lattice constant, 
(a) Minimum-energy structure: octahedral rotations (black) 
oftypeRno (*) and Rooi (+); polarizations (red) oftypePooi 
(o) and Pno (x). (b) Total energy (black squares) and c/a 
ratio (red triangles). 



cell), as in Fig. HJa). For larger lattice parameters, we 
see a transition at a — 3.91 A to an orthorhombic Ima2 
structure (space group 46, 10-atom primitive cell) with 
nonzero Pno accompanied by an Rno rotation of ~2°, 
combining the distortions of Figs, [ljb) and (d). Since 
there is no group-subgroup relationship between these 
two phases, the transition must be first order. These 
results are in rather good agreement with the previous 
first-principles calculations of Refs. and (assuming 
mixed domain phases are not considered), even though 
those works did not include the octahedral rotations. 
This suggests that the small rotations that we predict on 
the higher-a side of the transition do not have a very pro- 
found effect on the transition in PTO. We also note that 
all of the first-principles results disagree with the zero- 
temperature prediction of a phcnomenological analysis^ 
where a window of r phase (Pooi & n d Pno both present) 
was predicted at T=0. 

As is well known, bulk PZO has a complicated Pbam 
40-atom unit cell built from a 2y/2 x \[2 x 2 enlarge- 
ment of the parent perovskite unit cell^ The structure 
is described by a combination of several distortions, the 
largest being an R-point oxygen octahedron rotation (of 
type a~a~c° in Glazer notation) and an antiferroelectric 
(n/2a, w/2a, 0) E-line mode^S The latter is not allowed 
within the supercell used in our calculations. However, 



in bulk PbZr03, the pattern of atomic displacements for 
this mode is dominated by antipolar displacements of the 
Pb atoms in the [001] plane, and it is unlikely to survive 
in our [1:1] supcrlattice, in which the Pb atoms are shared 
with the PZO layer. Therefore, since our immediate goal 
is to provide a reference for the PTO/PZO superlatticc 
calculations, we continue to restrict our calculations to 
structures realizable within our 20-atom simulation cell. 
After considering a wide variety of structures consisting 
of different combinations of polar and rotational distor- 
tions, we find that the monoclinic Alal phase (space 
group 9, Clcl in the conventional setting, 10 atoms/cell), 
is lowest in energy throughout the range of in-plane lat- 
tice constants considered, as shown in Fig. [3J All four 
of the distortions Pooi, Pno, R-ooi, and Rno shown in 
Figs. HJa-d) are present in this phase. 

B. PTO/PZO superlattice 

We now turn to the case of the PTO/PZO superlattice. 
Compared to the relatively simple behavior for pure PTO 
and pure PZO, the phase sequence for the PTO/PZO 
superlattices presented in Fig. 0] is quite rich. 

We first remind the reader of the previous results 
of Bungaro and Rabe^ who confined their considcra- 



5 



tions to structures that do not include octahedral ro- 
tations. Within this constraint, these authors found a 
sequence of two phase transitions, first from a tetragonal 
P4mm structure ( "c phase" ) for compressive strains a < 
3.96 Ato a monoclinic Cm structure ("r phase") at inter- 
mediate strains, and finally to an orthorhombic Amm2 
structure ("aa phase") for tensile strains a > 4.12 A. Our 
results, presented in Fig. 21 show a qualitatively similar 
behavior for the polarization, although all of the phases 
we observe include octahedral rotations, and the phase 
boundaries (vertical dashed lines in Fig. @| have shifted 
to make the region of monoclinic phase wider. We give 
a detailed presentation of our results in the following, 
first in the region of compressive strain, and then in the 
tensile region, where the behavior is especially subtle. 

1. Compressive strain region 

Starting at the smaller lattice constants a < 3.890 A, 
we find a tetragonal P4bm structure (space group 100) in 
which Pool and Rooi are both present (20-atom primitive 
cell). The octahedral rotation angles for Rooi in the PTO 
and PZO planes are quite different, with the PTO rota- 
tions remaining rather small, ~1°. (Throughout the fig- 
ure, only the PZO rotations are plotted.) When the lat- 
tice constant exceeds 3.890 A, we find that Pn and Rno 
turn on simultaneously in a second-order transition, lead- 
ing to a monoclinic Plal structure (space group 7, Plcl 
in the conventional setting, 20-atom primitive cell). This 
structure is very similar to the Alal monoclinic structure 
obtained for PZO in Fig. [3l although the inequivalence 
of the PTO and PZO layers makes the Plal structure 
more complicated. The quoted value of a cr ;t=3. 890 A was 
obtained using the approach described in Sec. Ill CI but 
applied in the compressive-strain region, and using the 
distorted structure at a=3.925A to define the distortion 
eigenvector. 



2. Tensile strain region 

Next, as the lattice constant increases and approaches 
a value of about 4.19 A, we find that the amplitudes 
of Pool and Rooi diminish and approach zero together. 
Again using the method of Sec. Ill CI we find that these 
mode amplitudes disappear in a second-order transition 
at a=4.171 A, leading to an orthorhombic P2an phase 
(space group 30, Pnc2 in the conventional setting) in 
which only Pno and Rno are present (20-atom primitive 
cell). This structure is similar to the Ima2 phase seen 
at larger lattice constants for PTO in Fig. [31 except that 
the PTO and PZO planes are inequivalent. 

However, at a marginally larger lattice constant 
(a=4.193A), a first-order transition occurs to an or- 
thorhombic P2\am phase (space group 26, Pmc2\ in the 
conventional setting, 20-atom primitive cell). The value 
of a C rit=4.193 A was obtained from the crossing of curves 




FIG. 4: (Color online) Calculated properties of PTO/PZO 
superlattice as a function of epitaxially constrained in-plane 
lattice constant, (a) Minimum-energy structure: octahedral 
rotations (black) of type Rno (*) and Rooi (+); polarizations 
(red) of type Pooi (°) and Pno (x); antipolar distortions of 
type Aj 10 (filled red circles), (b) Total energy (black squares) 
and c/a ratio (red triangles). Heavy and light vertical dotted 
lines indicate phase boundaries determined from our work, 
and from Ref. 0, respectively. 



of computed energy vs. lattice constant as computed in 
the Plal and P2\am symmetries. The P2\am phase is 
one in which Rno disappears, Pno is approximately un- 
changed, and Aj 10 and Rooi distortions appear. A no , 
denoted by the solid red dots in Fig. |U consists of dis- 
placements (largest for Pb, then Ti, then Zr, but also 
involving O atoms) that are parallel and antiparallel to 
(110), patterned in such a way that no net polarization 
develops along this axis, as shown in Fig.flje). (Such an- 
tipolar distortions are also symmetry-allowed in the mon- 
oclinic Plal phase, but we find them to be very small, 
of order 0.02 A, suggesting that they are induced there 
by intermode couplings as will be discussed in Sec. lIVBl ) 
Contrary to the usual case, the signs of Rooi are the same 
in the PTO and PZO planes, with the rotations in both 
layers being small (~1°). The fact that the displace- 
ments associated with A no are much larger than those 
associated with R oi suggests that A no and Pno consti- 
tute the primary instabilities; when these modes freeze 
in, they induce an Rooi distortion, as will become clear 
from the symmetry analysis to be given in Sec. IIVI 

The fact that the second-order and first-order transi- 
tions occur at almost the same lattice constant has to be 
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regarded as a coincidence. The window of P2an phase 
found in our calculation, 4.171 A< a <4.193A, is in fact 
so narrow that we cannot be confident of its existence at 
all. Small changes in the technical details of the calcula- 
tions, such as the exchange-correlation potential or choice 
of pseudopotentials, could erase or widen this window of 
stability. The inclusion of c-axis tilting, as discussed in 
Sec. Ill Dt also tends to reduce or eliminate the window 
of Plan phase. We can only conclude that the Plan 
phase is competitive with the other two phases near the 
transition from the P4bm to the P2\am phase, and may 
appear in a narrow window between them. If the Plan 
phase is not in fact present, then there would be a direct 
transition from PAbm to P2\am. (Despite the presence 
of a group-subgroup relation, the transition would be first 
order because of the discontinuity in Rno-) 

IV. DISCUSSION 

This system illustrates several interesting aspects of 
complex perovskite oxides. Perovskite oxides generally 
are characterized by multiple unstable modes in the cu- 
bic structure, which couple and compete to produce low- 
energy states. The instabilities also each separately cou- 
ple to strain. Changing strain not only produces or elim- 
inates instabilities, but it can also change the balance of 
competition and produce first-order phase transitions be- 
tween distinct phases without a group-subgroup relation. 
Here, we see that combining two different perovskite con- 
stituents, each with characteristic unstable modes and 
strain dependence, can lead to a great increase in rich- 
ness of the observed phases. 



A. Parent phases 

In the present case, the computed phases in the two 
constituent materials are relatively simple. The domi- 
nant instability in cubic PTO is the zone-center polar 
mode, which couples to strain so that elongation is fa- 
vored along the polar direction. Compressive epitaxial 
strain thus favors Pooij while tensile epitaxial strain fa- 
vors Pno- In contrast to the proposed pressure-induced 
tctragonal-monoclinic-rhombohcdral phase sequence^ 
the absence of a monoclinic phase with rotating polariza- 
tion, or equivalently the first-order character of the tran- 
sition, can be understood as arising from the strong cou- 
pling of the polarization to tetragonal strain, as demon- 
strated by the quartic anisotropy of the polarization en- 
ergy given in Table III of Ref. [27|. Additional symmetry 
breaking for tensile strain results from the rotational in- 
stability Rno- While compressive strain is expected to 
favor the rotational distortion Rooi, no instability occurs 
in the range of strain considered here. 

In contrast, cubic PZO has multiple instabilities that 
occur together, yielding a monoclinic phase under the 
constraint of the superccll used for the first principles 



calculations. This phase is strongly favored and ap- 
pears for all values of strain considered, with compres- 
sive strain favoring the modes Pool and Rooi and ten- 
sile strain favoring the modes Pno and Rno, but not 
with such a strong dependence as to eliminate any of 
the modes at the highest strains considered. Also in 
contrast to PTO, the polarization freely rotates in the 
(llO) plane with changing epitaxial strain; this was al- 
ready apparent in the first-principles results for strained 
PZO neglecting rotations, which included an intermedi- 
ate monoclinic phase between the phases with normal 
and in-plane polarizationji^i 



B. PTO/PZO superlattice 

We now turn to the superlattice phase sequence. The 
point of closest resemblance to the phases of the con- 
stituents is the monoclinic phase for intermediate strain, 
which shows the rotating polarization and accompanying 
rotations observed in PZO. However, the strain depen- 
dence of the mode amplitudes is much stronger. With 
increasing compressive strain, Pno and Rno disappear, 
yielding a phase with only Pool an d Rooi- With increas- 
ing tensile strain Pool an d Rooi disappear, yielding a 
phase with only Pno and Rno, which at slightly larger 
in-plane lattice constant undergoes a first-order transi- 
tion to a phase dominated by Pno, Aj 10 and Rooi- 

A symmetry analysis of the Taylor expansion of the 
energy around a high-symmetry reference structure in 
the amplitudes of the relevant distortion modes is use- 
ful for understanding the mode contents of the various 
phases observed. For the pure compounds, the high- 
symmetry reference structure is taken as the tetrago- 
nal P4/mmm structure with 5 atoms per cell and all 
atoms at high-symmetry positions. For the superlattice, 
the high-symmetry reference structure is the 10-atom 
P4/mmm structure. We consider the distortion modes 
shown in Figure 1; for the pure compounds there is an 
extra translation symmetry relating the two B cations 
and the modes shown are either even (Pool; PiiOj A-iio) 
or odd (Rooi j Rno) under this translation. 

For pure PTO, the analysis is simple. For compressive 
strain, the Pool mode is unstable and is the only mode in- 
cluded in the lowest-energy phase. For tensile strain, the 
lowest-energy phase includes both Pno and Rno distor- 
tions. In previous worki we have seen that the Pno mode 
is unstable for tensile strain. Rno is thus nonzero either 
because it is independently unstable, as would be favored 
by tensile strain, or because the biquadratic coupling be- 
tween Pno and Rno is cooperative and large enough to 
destabilize Rno. 

For pure PZO, the four modes Pool, Pno, Rooiand 
Rnoare included in the lowest-energy phase. It is reason- 
able to assume that the polar distortions Pool and Pno 
are unstable throughout this range, with instability in 
Rooi promoted by compressive strain and the instability 
in Rno promoted by tensile strain. However, the rota- 
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tional modes need not be independently unstable. Bi- 
quadratic coupling to the polar modes could produce the 
effective instability. Furthermore, there is a fourth order 
invariant RnoPnoR-ooiPooi in the energy expansion that 
ensures that if three of the four distortions are present, 
the fourth will be induced through an improper mecha- 
nism even if it is not unstable. 

The effect of this fourth-order coupling is particularly 
striking in the superlattice, as it is responsible for the 
coupled appearance/disappearance of Pno and Rho at 
a=3.890 A and the coupled appearance/disappearance of 
Pooi and Rooi at et=4.17lA. The fact that these modes 
appear and disappear together might at first seem sur- 
prising, since in the high-symmetry reference structure, 
these two distortions have different symmetry and the in- 
troduction of one mode into the reference structure would 
not induce the other. However, with the lowered symme- 
try due to the presence of Pooi and Rooi, the two modes 
Pno and Rno belong to the same irreducible representa- 
tion and therefore are expected to occur together. This is 
formally equivalent to the presence of the fourth-order in- 
variant RnoPiioR-ooiPooi, discussed above. For example, 
in the compressive-strain phase with Pooi and Rooi, if an 
instability with respect to Pno or Rno develops, Rno or 
Pno will be linearly induced through this fourth-order 
term, producing the intermediate monoclinic phase and 
ensuring that the two modes appear/disappear together. 

It is in the extension of this symmetry analysis to in- 
clude the A no distortion that differences between the 
pure compounds and the superlattice arise. In the pure 
compounds, the Aj 10 distortion is a symmetry-breaking 
distortion that would appear only if it were unstable, 
which is not the case for the strain range considered. 
In the superlattice, there are two third-order couplings 
RooiPiioAi 10 and RiioPooiAjxo, not allowed in the pure 
compounds due to the additional translation symme- 
try. In the intermediate-strain monoclinic phase, both 
of these terms would linearly induce an Aj 10 distortion, 
which is indeed present at small amplitude in the com- 
puted structure for this phase. 

In contrast to its minor role in the intermediate-strain 
monoclinic phase, Aj 10 appears to be a strong indepen- 
dent instability in the superlattice at higher tensile strain. 
This can explain the first-order transition at a=4.193A. 
The approximate continuity of Pno through the first- 
order transition suggests that Pno is by far the domi- 
nant instability at higher tensile strain. Energy lower- 
ing from this polar phase can be produced in two ways: 
through freezing in an Rno instability to obtain the Pnc2 
phase, or by freezing in an Aj 10 instability, which induces 
a small Rooi distortion through the third-order invariant 
Aj 10 RooiPno and yields the Pmc2\ phase. The key to 
the first-order transition is the competition of these two 
instabilities in the presence of Pno- 



C. Growth of PTO/PZO and related superlattices 

We have attempted the growth of [1:1] PTO/PZO su- 
perlattices using pulsed laser deposition (PLD), but with- 
out notable success. Typically, reflection high-energy 
electron diffraction (RHEED) indicates that an initial 
layer-by-layer growth rapidly gives way to a 3D growth 
mode. Our best result to date is for growth on a NdSc03 
(NSO) substrate, where we observed five successful coher- 
ent repetitions of the [1:1] PTO/PZO unit before the 3D 
growth mode took over. These results were obtained us- 
ing the deposition settings that were optimized for PTO 
by Catalan et ali22 The settings for PTO were also used 
for the PZO growth. 

Previous work on SrTKVSrZrO^ and KTa0 3 - 
/KNbOg^i systems has shown that the synthesis of B- 
site superlattices is possible. We speculate that the de- 
velopment of new substrates (or new substrate treat- 
ments) might help extend the layer- by- layer growth 
of PTO/PZO, especially for larger lattice parameters 
(a=4.0-4.2 A). One option could be the new treatment 
developed for scandates (including NSO) by Kleibeuker 
et al^ 

Even if the extended growth of PTO/PZO superlat- 
tices remains elusive, this work may still provide useful 
guidelines for designing new ferroelectric materials based 
on B-site superlattices. This can be done either by creat- 
ing new low-symmetry phases by making use of fourth- 
order coupling terms such as Pooi Rooi PnoR-iiOj or by 
inducing improper fcrroclectricity via a third order term 
such as RooiPiioAi 10 or RnoPooiAno- For example, 
in the latter scenario, a system with an Aj 10 instabil- 
ity and an Rooi or Rno rotational instability should be 
an improper ferroelectric. For A-site superlattices, recent 
work on the PTO/SrTi0 3 system^ demonstrates that it 
is indeed possible to synthesize superlattices that induce 
improper ferroelectricity, even if the details of the sym- 
metry analysis are somewhat different in that case as is 
described by Rondinclli and FennieJ^ 

V. CONCLUSIONS 

In summary, we have presented the results of first- 
principles calculations of the structural phase transitions 
in [1:1] PTO/PZO superlattices, as well as in the parent 
compounds, as a function of in-plane epitaxial strain. A 
symmetry analysis was used to clarify the kinds of struc- 
tural distortions occurring in each phase and to illumi- 
nate the nature of the transitions between phases. 

While the PZO system remains in a single monoclinic 
phase over the whole region of lattice constants studied, 
and the PTO system has a single transition from a tetrag- 
onal to an orthorhombic phase, we find a significantly 
more complex phase sequence for the PTO/PZO super- 
lattice. Specifically, with increasing lattice constant, we 
find a sequence of three transitions, first from a tetrag- 
onal to a monoclinic phase, then to a first orthorhom- 
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bic phase, and finally to a second orthorhombic phase. 
The first orthorhombic phase occurs over such a narrow 
range of in-plane lattice constant that we cannot be con- 
fident of its existence in real PTO/PZO superlattices. 
The first three of these phases all contain strong polar 
and octahedral-rotation distortions, while the fourth has 
strong in-plane polar order, substantial antipolar distor- 
tions, and weak rotations. 

Our results indicate that polarizations and rotational 
distortions are strongly coupled in the PTO/PZO super- 
lattice. For example, for the monoclinic phase, in which 
four kinds of mode distortions are present simultaneously, 
we find these modes to be coupled such that the presence 
of any three of them will automatically induce the fourth, 
giving rise to a surprisingly rich behavior. Thus, we find 
that the inclusion of the octahedral rotations is crucial 
for an accurate determination of the phases and phase 
boundaries. 

While the synthesis of short-period perovskite super- 
lattices has the potential to lead to a wealth of new func- 



tional materials with improved properties, a thorough 
understanding of the interplay between epitaxial strain 
and spontaneous structural distortions will be needed to 
guide the search for such materials. The present work is 
a step along the way to this goal, and we hope that it 
might help identify new materials with enhanced ferro- 
electric and piezoelectric properties for next-generation 
applications. 
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